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Abstract 

Number-conserving cellular automata (NCCA) are particularly interesting, both 
because of their natural appearance as models of real systems, and because of the 
strong restrictions that number-conservation implies. Here we extend the definition 
of the property to include cellular automata with any set of states in Z, and show 
that they can be always extended to "usual" NCCA with contiguous states. We 
show a way to simulate any one dimensional CA through a one dimensional NCCA, 
proving the existence of intrinsically universal NCCA. Finally, we give an algorithm 
to decide, given a CA, if its states can be labeled with integers to produce a NCCA, 
and to find this relabeling if the answer is positive. 
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1 Introduction 



A cellular automaton (CA) is a discrete dynamical system, where the nodes 
("cells") of some regular lattice (usually Z d ) are mapped to a finite set of 
states; the state of a cell at a time t + 1 is determined by a local function 
that takes at inputs the states of the cell and its neighbors at t. Cellular 
automata have been widely used as models of dynamical systems in which the 
behavior is determined by local interaction between spatially fixed elements. 
An interesting particular class of CA is the class of number- conserving CA 
(NCCA): roughly speaking, these are CA where the states are represented as 
numbers, and the sum of the states over all cells remains constant when the 
states are updated. This conservation may help to prove some properties of 
the dynamics, and can be usually related to the conservation of some quantity 
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in the system that is being modeled. In CA traffic models ([13]), for instance, 
states are interpreted as the number of indestructible particles located in a 
cell. In fact, an interpretation in terms of particles can be given to any NCCA 
[10]. 

In [3], Boccara and Fuks give a necessary and sufficient condition for a one- 
dimensional CA of two states to be number-conserving, and study all the 
NCCA with neighborhoods {/,..., r}, I + r < 4. In [4], they give a necessary 
and sufficient condition that holds for one-dimensional CA of any number of 
states, and use it to study all the three-state NCCA for l + r < 2. In [7] Durand 
et al. formalize three different definitions of number-conservation and show 
their equivalence. They write the generalization of Boccara's condition to two 
dimensions, and hint on the (i-dimensional case; the decidability of number- 
conservation is thus proved (provided that the numeric value of the states is 
given). Durand et al. also give some examples of NCCA in several classes of 
one-dimensional CA, intersecting the classifications of Kurka ([9]) and Braga 
et al. ([5,6]), and prove the emptiness of the remaining classes. In [11] Morita 
and Imai prove the universality of the class of number-conserving reversible 
partitioned CA; in [12] Morita et al. embed a simple general computer in 
a reversible, number-conserving two-dimensional partitioned CA. We must 
remark that, despite the title of [11], these articles do not settle the universality 
of NCCA; partitioned CA can be recoded as normal CA, but the recoding does 
not, in general, preserve number-conservation. 

In section 2 we give the definition of NCCA and recall the necessary and suffi- 
cient conditions of [4] and [7]; we generalize the definitions and the conditions 
to allow any finite subset of Z as the set of states. In section 3, we show how any 
NCCA with such set of states can be extended to a NCCA with a set of states 
of the usual form {0, 1, ... , q}. In section 4 we show that any one dimensional 
CA may be simulated by a one dimensional CA. This implies that NCCA are 
capable of universal computation, and the particular form of the simulation 
implies the even stronger property of intrinsical universality. Finally, section 
5 addresses the question of deciding, for any CA, whether its states can be 
relabeled with different values in Z to make the CA number-conserving, thus 
showing the decidability of the number-conservation property in a wide sense. 



2 Definitions, previous results, and generalization 

The definition and results for NCCA have so far assumed a set of states of 
the form S — {0,1, ... , q}. This makes sense for some applications, but in the 
general case, there is no reason for restricting the class this way: we will define 
them for any finite S C Z. As we will see, this doesn't change the previous 
results, and turns out to be useful in the later sections. 
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Cellular automata: A cellular automaton F in the d- dimensional space is 
formally described by a tuple F = (d, S, N, f), where d G N is the dimension, 
S is a finite set of states, N C Z is a finite neighborhood, and / : S N — > 5 
is a /oca/ transition rule 1 . Since the neighborhood can always be enlarged by 
ignoring additional neighbors, we may assume it to be an hypercube, say, 

N = {(a 1 ,...,a d )dZ d :-l i <a i <r l V*} 

for some non-negative integers l±, r±, . . . , Id, r&. A configuration is an element 
c G S z . The dynamics of the system is discrete in time, and at each time step, 
the current configuration c* G S zd determines the next one, c t+1 , through 

tfti,...,i d ) = f( ct \(h,...,i d )+N) 

In this way the local function induces a global one, which we will denote with 
the name of the automaton: F : S zd — > S z ° '' . 

For the next definitions we follow Ollinger [14]; his formalization of intrinsic 
universality follows an idea presented in Albert and Culik II [1] and in Bartlett 
and Garzon [2]. 

Sub-automaton: Let F and G be two cellular automata with sets of states 
Sf and Sq respectively; we say that F is a sub-automaton of G if there is an 
injective map (f> : Sf — > Sq such that (f> o F(c) = G o 0(c) for all c G S F . In 
other words, F is just the restriction of G to some configurations (up to state 
relabeling). 

Scaling: Here we consider the case with d = 1. Let a be the shift function, 
i.e., the function a : S" 1 — > S z such that cr(c)j = For some strictly positive 
integer m, the packing map o m is the function o m : S z — > (S" 71 ) 2 such that 
o m (c)j = (c m j, . . . jCmi+m-i). Notice that both functions are bijective. If F is 
a CA with states S, a (m, n, fc)-rescaling of F is a cellular automaton F ( - m ' n ' k ' ) 
with states S m which verifies F< m - n ' fc >(c) = o- fe o o m o F n o (o m )" 1 (c) for all 
c G (S" 1 ) 2 

Simulation and universality: We say that a CA F simulates a CA G if 
G is a sub-automaton of some rescaling of F. A CA is said to be intrinsically 
universal if it simulates any other CA. This definition, thanks to the restricted 
definition of simulation, is stronger than the usual definition of computational 
universality, which asks for the ability to simulate a universal Turing machine. 

1 We restrict the definitions to what we use; for general theory of CA, consult [8]. 
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Period: A period p G N d for a configuration c £ S z is a vector such that 
c i+p(c) = (k, Vi G Z '. it can be easily checked that a period is preserved trough 
the iterations of a CA: for any F = (d, S, N, /), a period p for c will be a 
period for c' = F(c). The expression < k < p will denote the set of vectors 
{0,...,pi} x ... x {0,...,p d }. 



Number conservation: Let Cp(d,S) be the set of all the configurations 
in S z that admit a period (the -spatially- periodic configurations); for each 
c G Cp(d,S) choose a period p(c). We say that a CA F = (d,S,N,f) is 
number- conserving iff 

0<fc<p(c) 0<fc<p(c) 



Durand et al. consider CAs with S of the form {0, ...,q}, for which they 
discuss three different definitions of number-conservation, and show that the 
three are equivalent. The first is the one we just gave; the second asks that 
the sum over all Z d be conserved, for all finite configurations (configurations 
where q = for all but a finite number of i G Z '). The third definition asks, 
for all ce S zd , that 



lim TiTTi = 1 where fj, n (c) = V 



Ci 



Necessary and sufficient conditions: 



Two previous results that we will use are the necessary and sufficient condi- 
tions for a CA to be number-conserving in one dimension (proved in [4]) and 
in two dimensions (proved in [7]). Both assume a set of states of the form 
S = {0, ...,q}, and explicitly include in the equation. In fact, the only 
property of which is used is its quiescence: we say that a state s is quiescent 
iff f({s} N ) = s. It follows directly from the definition of number- conservation 
that all the states of a number-conserving CA are quiescent. 

In order to show that G S is not required in those results, we will deduce 
the result for d = 1 again, without that condition. Consider F = (1, S, N, f), 
with iV = {— /, . . . , r}. Let n be n — I + r + 1, let a G S be any state, and 
take any (xi, . . . ,x n ) G S n . Consider the configuration c consisting of infinite 
repetitions of 

^lj *^2; ■ • • i ■Eni Oii ■ ■ ■ i CL 
. ' 

n-1 
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If we apply equation (1) to this configuration (for period In — 1), we obtain 



n-1 



(n-l)a + J2 x k = J2 /( a »- ,a,x u ... ,x k ) + ,x n ,a,. 



k=i 



k=i 



n—k 



k=l 



fe-1 



Replacing x\ = a in this equation, we get 



n-1 



77, 



a + J2 x k = J2 f( a i - j a > ^2, ... , a; fc ) + /(a, x 2 , ••• , x n ) 



fe=2 



fe=l 



n-fc+1 



n 



a 



fc=2 
n-2 



fe-1 



^ /(a, ... , a, x 2 , •■• , Xfc+i) + /(a, z 2 , ... , x„) 



fe=i 



n—k 



+ J2f(x k ,- » I »» a '- ' a ) + 



fc=2 



fe-1 



where we used the fact that /(a, . . . , a) = a. Taking the difference between 
the two last equations, we have 



n-2 



xi - a = /(flv^a^i, ••• ,x k ) - f(a 1 ^a,x 2 , - ,x k+ i) 

^ ^ n— fc n— fe 

+f(a, xi, ... , x n _i) - f(a, x 2 , ... , x n ) + /(xi, ... , x n ) - a 



i.e., 



n-1 



/(xi, ... , x n ) = x 1 + fia^-jQ, x 2, - , a:fc+i) - /(fV-^a, a?i, - , a*) (2) 
fc=i 



n—k 



n—k 



This is a necessary and sufficient condition. If F verifies equation (2) for all 
(xi, . . . ,x n ) G S n , then it verifies (1): the terms in the brackets will cancel 
when the sum ranges over the whole configuration, leaving just the sum of 
states before (on the left) and after (on the right) the application of F. Hence, 
the theorem reads: 

Theorem 1 Let F = (l,S,N,f) be a CA with S C Z, N = {-/,..., r}, 

n — I + r + 1. Let a be any state in S . Then, F is number-conserving if and 
only if f verifies equation (2) for all (xi, . . . , x n ) G S n . 
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The analogous condition for two dimensions was proved in [7] (for a = 0). The 
explicit writing of the proof is cumbersome, but the idea is the same as before: 
this time we consider a configuration with a matrix ( 2^)1=1,.. .,^=1,...^ (the size 
of the neighborhood) surrounded by a's, and write the necessary condition. 
Then, instead of subtracting the evaluation with x\ = a as before, we subtract 
the evaluation with x± t , = a and the evaluation with x, t i = a, and add the 
evaluation with both = a and x %t i = a. It is easy to see how the procedure 
can be further modified to get necessary and sufficient conditions for d > 2. 
For two dimensions, the resulting condition can be written as follows: 

Theorem 2 Let F = (2, S, N, f) be a CA with S C Z, N = {-l u . . . , rj x 
{—I2, • • • , r 2 }, n — l± + r\ + 1, m — l 2 + r 2 + 1. Let a be any state in S. Then 
F is number- conserving if and only if, for all ...x m , n ) G S nm , it satisfies 

f(M 1Xm%n ) = 

m—l n—l 

xi,i + £ /(M 2 , M+1 , n ) - /(M liM , n ) + £ /(M 1)2 , mJ+1 ) - f(M 1M ) 
i=i j=i 

m—l n—l 

+ E E f( M i,2,i,j+i) + f(M 2 ,i, i+ i,j) - /(Mi )MJ ) - f(M 2t 2,i+i,j+i) (3) 
i=i 3=1 

where Mt,l,b,r represents a matrix filled with a 's but for the bottom left corner, 
which is occupied with (xij)i=T,...,B,j=L,...,R- 

3 Extension from S C Z to {0, . . . , q} 

The theorem of this section shows that the introduction of general S C Z as 
possible sets of states does not change the class of NCCA in any dramatic way: 
anything that can be seen in a NCCA with S e Z, can be seen in a NCCA 
with S of the form {0, ...,q}, with the appropriate initial conditions. The 
theorem is stated and proved in one dimension, but we remark that versions 
of it for higher dimensions can be proved with similar arguments. 

Theorem 3 Let F be a NCCA with states S C Z and neighborhood of size 
n. Then F is a sub-automaton of a NCCA F with states S — {0, . . . , max S — 
minS"} and neighborhood of size In. 

Proof. Let F be a NCCA, F = (1, S, N, /), with N = {-/, . . . , r}. By sub- 
tracting minS" from all the states in S, we can assume that S C {0, . . . , M}, 
with M = max S — min S and G S, M G S. In addition, without loss of gen- 
erality, we may assume / = 0: if I > 0, we can apply the result to F' = Foa~ l , 
obtaining F', and then take F = F' o a 1 . 
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We first note that for any c G S z , 

i 

]T f( Cj ,...,c j+r ) < rM + Ci (4) 



j=i-r 



To prove this assertion, notice that for this to be false, c\ must be strictly 
smaller than M; but in that case, we can change c, putting q = M, and 
we know, from the number-conservation, that the sum in F(c) must increase 
accordingly. Since the only cells that can notice the change are those that 
"see" the state of i, we have that Y?j=i- r F{ c )j m ust increase in M — C{. Since 
the new sum is bounded by (r + 1)M, we get (4). 

The procedure to get F(c) from a configuration c is the following: 

C; for Cj G S 



~ 



c" for a G S 



(1) Viez,4 

(2) c" = F(c') 

(3) Vi G Z , c° 

" < + c, 

(4) For k— 1 , . . . , r 

Vi G Z , = max{0, c? -1 - M} 
VieZ,c* = c?- 1 -e? + e*. 1 

(5) F(c) = c r 



So, we first remove q from each "invalid" position i, obtaining d G S" 2 . We 
apply F (which is a number-conserving operation), and then we put q back 
in its position, recovering the original sum of the states). At that point, some 
c° may be greater than M. To correct this, the surplus at each site is pushed 
to the left, and this is done r times. 

We claim that c r verifies < c\ < M, for all i. Think of the q that are added 
in step 3 as particles labeled with i, and assume that each time an surplus is 
moved to the left (in step 4), preference is given to the particles with a higher 
label. Equation (4) assures that in the r sites to the left of i, there is enough 
place to accomodate the q particles: 

( r + l) M - £ c" = M + rM- ]T ./(/ > M > Cl 



(remember that c- = for the i we are considering). Notice that no particle 
is moving more than r steps. Therefore, the final state c\ is determined by 
c'[ +k and c i+ k for k — 0, . . . , r; i.e., it is determined by c i+ k for A; = 0, ... , 2r. 
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To define /, we just consider each (xo, ■ ■ ■ , x 2r ) G S 2r+1 , apply the preceding 
procedure to the configuration c defined by Cj = x« for i = 0, . . . , 2r and 
q = elsewhere, and set f(x , . . . ,x 2r ) = c 7 . If (x , . . . ,x 2r ) £ S 2r+1 , then 
d = c in the procedure, there are no surplusses to move, and f(xo, . . . , x 2r ) = 
f(x ,...,x r ). □ 



4 Universality 

Theorem 4 Let F be a CA with q states and neighborhood of size n. Then F 
can be simulated by a NCCA G with states {0, . . . , 2q + 1} and neighborhood 
of size In. 

Proof. Let F be F = (1, S, N, /), with 5 = {1, . . . , q} and TV = {-I, . . . , r}. 
We define F — (1, <S, N, /) with 

5 = {-g, . . . , -1, 0, 1, . . . , q} , AT = {-21 - 1, . . . , 2r + 1} 

and f(a-2i-i, a 2r +i) given by 

I/(a_2i, a-2«+2, ■■• , a 2r ) if a 2i = -a 2i+ i > , -/ < i < r 

f(-a-2i, -GL2H-2, ■■• , ~a 2r ) if a 2i = -a 2 i-i < , -/ < i < r 
ao otherwise 

The idea is the following: each cell is split in two, and the new cells are occupied 
with a negative and a positive copy of the state the cell had. The iterations 
will follow the original CA rule, both on the negative and on the positive cells, 
and the sum will remain constant (zero). 

With the definition given above, a positive cell will change its state only 
if it sees a "correct configuration" around it (a configuration of the form 
a, —a, b, —b, . . .), and sees that its right neighbor is seeing it too; similarly, a 
negative cell will change only if it sees that its left -positive- neighbor will 
change in the same way (with the opposite sign). Hence, the only changes 
in the configuration are done in pairs of cells, and on each of this pairs the 
sum is conserved (and is 0). The CA F is then number-conserving. Through 
the injection a — > (a, —a) F becomes a sub-automaton of F^ 2 ' 1 ^ . To avoid 
"negative" states, we add q to all the states of F. □ 

Corollary 5 There are intrinsically universal one- dimensional NCCA. 
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Proof. The relation of simulation is a preorder [14], and there exist intrinsi- 
cally universal one-dimensional CA. □ 



5 Decidability 



Equations (2) and (3) show necessary and sufficient conditions for a CA in one 
and two dimensions, respectively, to be number-conserving; equations similar 
to them for higher dimensions are hard to write, but not to obtain. The de- 
cidabily of the property of number- conservation is thus proved, but with one 
important restriction: the numeric values of the states are taken as given. This 
is not a complete answer, since the labeling of the states in a CA is, in princi- 
ple, arbitrary. If a model produces a CA with a set of states defined as colors, 
letters, or even numbers, can we relabel the states with numbers such that the 
CA is number- conserving? If the answer is yes, it would be interesting, since 
the conserved quantity could be traced back to the original system, or could 
help to prove some property of the CA. In particular, in [10] it is shown that 
NCCA can always be interpreted in terms of the interactions of indestructible 
particles; it would be interesting, for a system with states, say, blue, yellow 
and red, to know if the dynamics can be expressed in such terms. 

If the set of states is required to be of the form {0, . . . , q — 1}, where q = 
\S\, then the answer to the problem is easy: just consider the q\ possible 
permutations of the labels, and check, for each of them, if the CA is number- 
conserving (using equations (2) or (3)). But this requirement is arbitrary: 
perhaps what we are seeing can be interpreted in terms of particles, but there 
is some quantity of particles which happens to never occur in a cell, and is 
therefore absent from our current set of states. In general, we would like to 
know if we can relabel the states of the CA with l^l different elements of Z, 
so as to make the CA number- conserving. 

Example 1 Consider F = (1, S, N, f) with S = {a, b, c}, I — r — 3, and / 
defined by: 



f(x ,xi,x 2 ,x 3 ,X4,x 5 ,x e ) 



if 




= Xl 


= x 3 


= X4 


= a A X2 


= b) 


V 


(xi 


= x 2 


= X4 


= x 5 


= a A X3 


= b) 


V 


(X2 


= x 3 


= x 5 


= Xq 


= a A X4 


= b) 



X3 



It can be checked, using equation (2), that all the bijections <fi : {a,b,c} — > 
{0, 1,2} produce rules that are not number-conserving. But F does become 
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number-conserving if we relabel its states with 
0(a) = 4>{b) = 3 0(c) = 1 

It may be interpreted in terms of particles: if three particles are in a cell and 
the two neighboring cells in each direction are empty, then they separate; one 
goes to the left, one to the right, and one stays in the cell. 

Theorem 6 Let F = (d,S,N,f) be a CA. Then there is an algorithm to 
decide if there exists a relabeling that makes F number- conserving, and to find 
it if the answer is positive. 



Proof. Write S as S = {s , si, . . . , s\s\-i}- We are looking for an injective 
function <f> : S — > Z, such that F, redefined with states (f>(S), becomes number- 
conserving. To find it solution, we will use the equation that gives the necessary 
and sufficient condition for the CA to be number-conserving; it will be equation 
(2), (3), or a form for higher dimensions, depending on the dimension in which 
F is defined. In [4] , the equation was used to list the possible NCCAs, taking 
the values of the function / as the unknown variables. This time, the unknown 
variables of the equation are not the evaluations of the rule, but the numeric 
value of the states: the necessary and sufficient condition to make the relabeled 
version of / number-conserving is, in the one-dimensional case (the other cases 
are analogous), that 

n-1 

(j>(f(x 1: ... ,X n )) = <j>(xi) + ^(/ (Scb ■■■ , S , X 2 , ... , S w -fc + i)) 

k ^ k 

-(f>(f(s , ... , S , Xi, ... , X n - k )) (5) 
k 

for all (x±,X2, ■ ■ ■ ,x n ) G S n , where n = \N\. (Notice that from theorem 1, 
the choice of so is arbitrary.) This is an homogeneous linear system of \S\ n 
equations and \S\ variables. The set of solutions will be a linear subspace 

V G IR' 5 '. If V — {0}, the algorithm gives a negative answer. If V ^ {0}, 
we must still find out if it contains solutions which have all the coordinates 
different (if not, then there is no injection). In other words, we must check if 

V \ (jjj^k Ejkj 7^ 0, where Ejk is the hyperplane Xj = Xk- But 

V\(\jE jk )=Q <=► v = vn(\jE jk )=\JvnE jk 

j^k j^k jy^k 

3j, k : V = V n E jk 3j,k:VC E jk 



where we use the fact that each V H E jk is a subspace of V, and a linear space 
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cannot be a finite union of proper subspaces [15]. The last condition can be 
easily checked by the algorithm, by adding the equation 4>{xj) = 4>{ x k) to the 
equation system and seeing if the space of solutions is the same. 

If V \ (jJj^k Ejk) 7^ 0! then there is a solution to the problem, and it can 
be found in finite time. We can now relabel A with 0, making it number- 
conserving; if a set of states of the form {0, . . . , q} is desired, theorem 3 may 
be applied. □ 

Note that the algorithm in theorem 6 is only needed when \S\ > 2. For l^l =2, 
if the CA can be made number-conserving, then it will be number-conserving 
for the two possible choices of <fi : S — > {0, 1}: if S = {a, b} C Z is a solution, 
then = {0, 1} is one too. 



6 Conclusions 

In this paper we develop several notions and results related to the class of 
number-conserving cellular automata, generalizing the definition and showing 
both the computational universality of the class and the decidability of the 
property that defines it (even for CA that are not initially expressed with 
numerical states). 

The generalization of the definition allows the CA to have any set of states 
in Z; we show that the necessary and sufficient conditions previously given 
for number-conservation do not require the previously assumed set of states 
of the form {0, 1, . . . , q}. Furthermore, the generalization does not change the 
class in any profound way: in section 3 we show that for a NCCA with set 
of states S C Z, the set of states can be completed to a contigous set of 
states, extending the rule while keeping the number-conservation property. 
Hence, any "generalized" NCCA can be seen as a "usual" NCCA, restricted 
to a subset of its configurations. The extension we show in theorem 3 requires 
an enlargement of the neighborhood; since we do not know examples where 
this enlargement is really needed, it may be the case that the result can be 
improved to an extension that preserves the neighborhood of the CA. 

Section 4 shows that any one-dimensional CA can be simulated by a one- 
dimensional NCCA; the construction is very simple and proves that members 
of the latter of the latter class exhibit intrinsical universality, a property that 
implies computational universality. 

Section 5 dealed with the question of deciding, for a given CA, whether its 
states can be labeled with integers, making it number-conserving. This ques- 
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tion turns out to be decidable, and we show how to find the relabeling, if there 
is one. The particular interest of this result is that it may help to reveal some 
conservative dynamics in an otherwise not-conserving CA; such a conservation 
may afterwards be interpreted in terms of indestructible particles [10]. 
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